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Most bright sources of the radiation in hard x-ray and gamma-ray regions are undulator sources 
and Compton based ones. These sources are ultimate for production of polarized positrons necessary 
for future linear colliders ILC, CLIC. We developed a novel method for evaluating the energy spec¬ 
trum of electrons emitting the undulator- and the inverse Compton radiation. The method based 
on Poisson weighted superposition of electron states is applicable for whole range of the emission 
intensity per electron pass through the driving force, from much less than unity emitted photons 
(Compton sources) to many photons emitted (undulators), and for any energy of the photons. The 
method allows for account contributions in the energy spread both from the Poisson statistics and 
diffusion due to recoils. The theoretical results were confirmed by simulations. The electron energy 
spectrum was used for evaluation of the on-axis density of photons and their coherency making use 
of the ‘carrier-envelope’ presentation for the emitting photons. The evaluated maximum coherency 
degree of single-electron radiation is evaluated to be proportional to the undulator spatial period 
and inversely to the energy of electrons, the number of coherently emitted undulator periods almost 
independent of the undulator deflection parameter. The results of our study are applicable both for 
the classical limit of classical undulator and for the quantum limit of Compton gamma-ray sources. 

PACS numbers: 41.60.Ap 


I. INTRODUCTION 

Most bright sources of the radiation in hard x-ray and 
gamma-ray regions are undulator sources and Compton 
based ones. These sources are ultimate for production of 
polarized positrons necessary for future linear colliders 

Both of the processes - undulator radiation and Comp¬ 
ton inverse radiation - exploit the same physical principle 
(radiation emission by electrons moving along sinusoidal 
or helical trajectory), and are known and theoretically 
described a few decades ago, see, e.g. M- 

Classical description of the undulator radiation pre¬ 
dicts maximal brightness attained on the undulator axis 
proportional to squared number of the periods of undula¬ 
tor magnetic structure, oc Ay due to spatial interference 
of the waves emitting by the electron. 

With increasing the electron energy and/or shorten¬ 
ing the undulator periods, the energy of undulator radi¬ 
ation advances to hard x-ray or gamma-ray region where 
the quantum effects in radiation emission will be ex¬ 
pected. The undulator radiation would resemble the in¬ 
verse Compton radiation, brightness is proportional to 
number of the laser photons with which the electron in¬ 
teracts, oc N \ as (no interference of individual photons is 
expected). 

In Q there was emphasized that under the quantum 
approach, instead of number of the undulator periods one 
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should account for number of photons emitted, which ac¬ 
tually much smaller than the number of periods. Nev¬ 
ertheless, since number of photons emitted directly pro¬ 
portional to number of periods, the result of [(J yields 
the same quadratic dependence of brightness as in the 
classical approach. 

In other words, these two approaches, the quantum one 
based on the Klein -Nishina electron-photon interactions, 
see [7||, and the classical approach yield different spectral- 
angular distribution of radiation. The quantum approach 
implies statistically independent quanta emitted, while 
under the classical approach all the quanta considered 
interfering with each other. 

We studied limitation of coherency of undulator radia¬ 
tion caused by recoils undergone by ultrarelativistic elec¬ 
trons emitting the photons. Due to the recoils, the energy 
of the electron is degraded and acquired a finite spread 
(due to a statistical nature of emission of the individual 
quanta). Under the classical approach, the undulator ra¬ 
diation is suggested either to keep the electron energy 
constant [1, (tj, or to adiabatically decrease (‘chirping ef¬ 
fect’ that can be eliminated by means of tapering the 
undulator, see 0)- 

Papers on quantum approach in undulator radiation 
predicted different energy distributions: in 0 there was 
obtained a brush like distribution depending on the num¬ 
ber of the undulator periods, in p the distribution is 
continuous depending on number of quanta radiated out. 
In a recently published paper [lj|] , the energy spectrum of 
electrons moving along an undulator is evaluated within 
the diffusion approximation. The energy distribution is 
continuous, some fraction of the initial delta ensemble 
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got acceleration (the authors commented it as somewhat 
nonphysical). A case of multiple Compton backscattering 
considered in 0 shows diffusive behavior of the electron 
spectra. 

The present paper is intended to cover the gap be¬ 
tween a classical undulator radiation of relatively low- 
energy electrons and a quantum one of high-energy elec¬ 
trons and the inverse Compton radiation (characterized 
by small number of quanta emitted per the pass through 
laser pulse). 

The paper is structured as follows. In the first sec¬ 
tion there described is a method to compute the energy 
spectrum of electrons passing along the undulator axis. 
The method is applicable for any number and any am¬ 
plitude of recoils. The second section comprises results 
of study on coherency of photons emitted along the axis 
by electrons with statistical parameters evaluated in the 
first section. A method of ‘carrier-envelope’ presentation 
is employed, which allows for estimation of the degree of 
coherency. A criterion for coherency degree is proposed. 
The third concluding section summarises the results and 
discusses possibility to attain coherency and thus high 
density of radiation in gamma-ray sources. 

Numerical examples explore the ILC baseline source 
of polarized positrons: energy of electrons 150 GeV, a 
helical undulator with the period 11.5 mm, active length 
127 m and the undulator parameter K = 0.92 (see El), if 
other values are not indicated. 



No. photons 


FIG. 1. Expected number of emitted photons for £ = 
0.1,1,10,50. 


length: 


fo(0 = e 5 • 


A vector of photon number probabilities (weights) can 
be defined with the components given by Eq. ©: 


f = {fn(0}- ( 2 ) 


II. ELECTRON SPECTRUM ALONG THE 
UNDULATOR 

A. General method 

Consider a certain pass length along the undulator 
(from the entrance to the coordinate z ), in which an elec¬ 
tron in average emits £(z) photons: 


B. Spectrum of electron emitted n photons 

Energy spectrum of an electron passing along an un¬ 
dulator can be presented as composed from the partial 
spectra corresponding to recoils caused by emission of n 
photons. We attribute the n-th electron state as average 
energy spectrum, normalized to unity, due to emission of 
n photons. Then the vector of states (in Hilbert space) 
can be constituted as: 


_ average energy loss 
average energy of photon 

The average number of emitted photons is increasing lin¬ 
early with 2 from £( 0 ) = 0 up to £(L„) = £* at the 
undulator exit. 

Number of emitted photons is supposed to obey the 
Poisson distribution law as being the statistically inde¬ 
pendent rare events (see e.g. |14l4la |l: 


fn(0 = r ; , n 

T(n + 1 ) 


( 1 ) 


where n = 1 , 2 ,... is the number of photons, T(n) the 
gamma function. 

Figured] presents probability of emitting n photons by 
the electron at different average number of photons £. 

Eq. |D yields that the fraction of particles has not 
emitted a photon exponentially decreased with pass 


F = {F n (r^o,K 1 X u )} (3) 


with normalization 


Ek (7) dy = 1 , 

where 70 is the initial energy (Lorentz factor) of the elec¬ 
tron at undulator’s entrance, K , A u the undulator param¬ 
eter and period, respectively. 

Emission of a photon with energy hio decreases the 
energy of electron: 

. frio 

7=7- 2 = 7 - 7 P h • 4) 

mc z 

It is induced spread in energy as well. 

Let W( 7 p h, 7 ) be a normalized to unity spectrum of 
emitted photons i.e. spectrum of recoils - (here and 
below we omit indication of the undulator parameters 
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A u ,iv). Then, in general, the components of the vector 
© can be sequentially calculated: 

F nh)= -Fn— 1(7 + 7ph)W / (7 P h,7 + 7 P h) dy p , (5) 

Mph 

starting with Fo ( 7 ) = 5(7 — 70 ) (5 is the Dirac delta- 
function) . 

Making use of centered the recoil spectra around the 
average value y p h = e + b with 6 ( 7 ) =< 7 P h >, Eq. © 
may be reduced to 

Fnil-b) = J F n _ 1 ('y + e)w(e,'y)de . ( 6 ) 

Thus n-th state is translation by b the convolution of 
71 —I state with the recoil spectrum. 

The aggregate spectrum of the electron having emitted 
in average £ photons is dot product of the state vector 
© and the vector of ‘weights ’ © 

■r(7,0 = = (7) 

ra =0 v ' 

For the undulator and Compton cases the recoil is 
small, the spectrum of recoil u>(e) is almost independent 
of 7 : b /7 « b /70 <C 1. For this case 6 ( 7 ) = 6 ( 70 ) is con¬ 
stant, the average energy of the n state and the variance 
are 

< 7 >ra = 7o -nb] 

M =< (7 + nb- 70) 2 >= 7 KJ 2 . 

Making use of these simplifications, the first two mo¬ 
ments of the aggregate spectrum read: 

< 7 >« = /«(£) J d 7 Fn ( 7)7 = 70 - M ; ( 8 ) 

<7 2 >i = J2fn(0 J d7F n (7)7 2 

= 7o + £(M - 2 7 o5) + £(1 + £)b 2 i (9) 

cr| =< 7 2 > — < 7 > 2 = £((X 2 + b 2 ) . ( 10 ) 

C. Estimations for undulators 

Let us consider, as an example, application of the de¬ 
veloped method to the case of ‘weak’ undulator, which 
is also equivalent to the Compton backscattering source. 
The spectral density of the undulator photons in the fun¬ 
damental harmonic is (see e.g. i): 

MC, K ) = G(K) [1 - 2C(1 - 0] e[l - C(1 + K 2 )} , (11) 

where G(K) is the normalizing factor, 0(x) represents 
the Heaviside theta function; ( is a reduced energy of 
photons; ( = e pho t /(2y 2 fiw u ) ■ 



FIG. 2. The first five analytical electron states, K 2 -C 1. 

For the case of ‘weak’ undulator, K 2 -C 1; G(K) = 
3/2, the centered spectrum of recoils reads: 

7 ) = g^ 3 ( 6 2 + e 2 ) 0(6 + e) 0 ( 6 -e) ( 12 ) 

where b = 7 qAc/A u with Ac being the Compton wave¬ 
length, A u the undulator period. 

Substituting CE2D into ©, we can successively derive 
the the electron states. The first five states are drawn in 
Fig. m 

The electron states possess following specific charac¬ 
teristics: 

• The average energy is < 7 > n = 70 — nb. 

• The base span is (< 7 > n — nb ) < 7 < (< 7 > n 
+nb ), no density above 70 and below 70 — 2nb. 

• The variance of n-state density is a 2 = na 2 with 
a 2 = lb 2 . 

• With increase of n the density shape gradually loses 
its individuality and tends to the Gaussian shape 
according to the Central Limit Theorem. 

It should be emphasized that the major contribution 
to the cumulative energy spread m comes from Pois¬ 
son superposition of the states £ 6 2 as compared with the 
width of the spectrum a 2 : 

M(0=£(M + M = £& 2 ( 1 + 2/5). (13) 

It also should be mentioned that for ID model consid¬ 
ered in El (suggested a spectrum of emitted photons be¬ 
ing of a delta-function shape, W{ 7 P , 7 ) = 5(y p — 7 ™ ax ) ) 
the spectrum of electrons derived by means of our method 
coincides with that obtained by authors of El : narrow 
picks distributed according to the Poisson law. 

Cumulative spectra of the first 5 states for relatively 
small £ are presented in Fig. [3] for large £ (all the states 
approximated by the Gaussian functions) in Fig. [I] 
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FIG. 3. Cumulative spectra for the different £’s. Superposi¬ 
tion of the first 5 states. 



FIG. 4. The central component, n = £, (higher), and the 
cumulative spectra (wider) for £ = 10, 30, 50; b = O.OOI 70 . 


D. Number of emitted photons 

Based on a brief description of the undulator radiation 
BUI, we estimate the length along an undulator axis 
to emit one photon in average by the electron. 

The total energy (in average) radiated by a single elec¬ 
tron traversing the iV u -period undulator is given by 

Pt = 7^ahN u ui u 'y 1 K 2 , (14) 

where N u is the number of undulator periods, w u = 
2 -7rc/A u ) the undulator equivalent frequency, A u the un¬ 
dulator period, K the the deflection parameter, a the 
fine structure constant. We suggest a helical undulator, 
for a planar one the energy is half of m- 

The average energy of the undulator photon is almost 
independent of the undulator parameter K and equal to 
half of maximum energy in the fundamental harmonic at 
K —> 0, see [f|. 

Substituting the average recoil energy b into (THD . we 
can derive the number of undulator’s periods v\ to emit 
one photon by the electron: 


3 

1 47TQ!iT 2 
For the case of K = 1, we have 

^heiica! _ 137 / 4 ~ 17 periods . 


(15) 


Accordingly, for the planar undulator zzf )lanar ~ 137/3 « 
46 periods. This estimation is in agreement with 9j 
where ~ 1/a periods. 

Transition to the planar undulator corresponds to re¬ 
placement K —» K / \[2 in m- For the Compton sources, 
the undulator spatial period should be substituted by 
A u —> Ai as /2 (for the head-on crossing of the electron 
with the laser pulse, see e.g. 0 )- 


In order to validate the analytical estimations, a 
Monte-Carlo code has been used to simulate the pro¬ 
cess of undulator radiation. The code, mainly intended 
for simulation performance of the Compton sources of 
hard radiation (see, e.g. jlTj), simulates the recoils of 
electrons due to a random emission of individual quanta 
with the spectrum m■ Fig. [5] represents spectra of elec¬ 
trons from simulations. Total number of particles 10000, 
the maximum number of quanta 800, the maximum en¬ 
ergy of photons e* = 0.002£’i n i t , the deflection parameter 
K <C 1. 

As it is seen from the figure, the recoil spectrum resem¬ 
ble the fundamental harmonic with additional low-energy 
tail of small intensity for £ <1. For the practical cases, 
£ 1 the energy spectra resemble the normal distribu¬ 

tion. 


III. FREQUENCY SPREAD AND COHERENCE 
A. On-axis density and coherence 

In the preceding section we regularized the electron 
kinetics: the energy spectrum of the electron passing 
along an undulator was derived as function of its po¬ 
sition, initial energy, and the undulator parameters: the 
field strength and the period. In the present section, we 
use two first moments of the spectrum, the average en¬ 
ergy and the spread (r.m.s. width of the spectrum), to 
evaluate a coherence degree of radiation emitting by the 
electron. The radiation will be presented as individ¬ 
ual pulses, uniformly distributed along the undulator, the 
distance between consecutive pulses is equal to T u /£l (^u 
is the undulator length). We will consider coherence of 
the photons in the fundamental harmonic emitting along 
the undulator axis - ID model. 
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FIG. 5. Simulated spectra of electrons. For £ = 0.1 - top left - the non-recoiled particles are not displayed. 


The degree of coherence will be estimated with the fol¬ 
lowing reasoning. The pulses are considered coherent if 
their field strengths are piled up, incoherent if the en¬ 
ergies of the pulses are summed. The density of energy 
provided by j pulses reads 


*>(j) 




dr , 


(16) 


If the pulses are normalized quadratically to unity, 

J ul(r)d t = 1 , 

then j pulses are considered fully coherent if power den¬ 
sity is proportional to j 2 , and fully incoherent if propor¬ 
tional to j. A degree of coherency may be attributed as 
(see, e.g. [13, [HI): 


m 


In D(j) 
In j 


(17) 


where u n {r) is the electrical field strength of n-th pulse, The degree of coherency m varies from T = 1 - 
r n is the offset. fully coherent pulses, energy density oc j 2 - to Tj = 0 
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fully incoherent, the maximum density oc j 1 . As it 
follows from the coherency is determined by the 

terms = / u m (r - r m )it n (r - r„) dr. Full coherency 
is equivalent to "frn,n = 1, incoherency to 7 m , n = 8 mn 
(Smn is the Kronecker delta). 


Since we suggested that the frequency spread in ra¬ 
diation is induced by the energy spread of the electrons, 
and the former resembles the normal distribution, we will 
consider gaussian pulses (packages). Introducing a nor¬ 
malizing factor P such that / u 2 (x, t ) d(x — ct) = 1, and 
substituting it into m, we get 


B. Regularization of emitted pulses 

The pulses generated by the electron passing a long 
undulator, L u /X u 1, may be represented with the 
“carrier-envelope” signal, see e.g. [22], where the car¬ 
rier frequency corresponds to the average energy of the 
electron, the pulse length is inversely proportional to the 
frequency spread induced by the spread of electron en¬ 
ergy. 

For a general form of a pulse, 


t ( T ; ko,cr) = 


a\pK [1 + exp(—fcgcr 2 )] 
_2 


exp 


-^2 ) cos ( fc 0 r) 


( 20 ) 


with t = x — ct and ko = 1/A carr i er . 

The interference term of two gaussian pulses, 


u(x,t ) 


\/27T J-c 


A(k) exp [—i(ui(k)t — kx)\ d k , (18) 


R(x; k, a) 


H exp 



cos (kx) , 


where x is scaled along the undulator axis, the photon 
pulse spectrum, A(k) is related to the pulse envelope with 
u(x, 0) via the Fourier transform (1191) . and vice versa. 

In a dispersion-free space, u(k) = kc, it reads: 


A{k) 



u(x, 0) exp[— ikx] da;. 


(19) 


H (<r, k ) 



At these conditions a pulse preserves its form and travels being the norm of a pulse, offset by r, possesses an exact 
at the speed of light c. form 

I 


7 m,n — / l?(x, <7 m )-/?(a; t, k ni o ji) da; — 


2cr„,cr n 


(^m + O-n) t 1 + ex P(-O-mfcm)] I 1 + exp(-(J 2 fc 2 )] 


T i.(km-k n r+T 
2(°’m - l -cr n ) 


COS 


km^m + k n al 


n (T~( fcm+fen) 

2 (v'tn+v'A) 


COS 


2 1 


k rn a m k n (J n 




( 21 ) 


In practical cases a photon pulse in undulators comprises many periods, ak < 7V U . Therefore exp(— a 2 k 2 ) 
exp(— N 2 ) 1 and the interference term (12T1) can be simplified: 



exp 


~ k n ) 2 + T 2 
2 (^m + <?%) 


COS 


kr„O rn 7 k n (T n 


( 22 ) 


Expression (1221) allows for computation the power den- C. Connection to the undulator radiation 

sity m and consecutively the coherence factor. The 
term m consists of the exponential factor depending of 
widths of pulses, and an oscillating cosine factor respon¬ 
sible for variation of the carrier frequency between two 
pulses. 

Within adopted approximation of small recoil, neglect¬ 
ing an initial energy spread in electrons and the frequency 
spread due to the maximal pulse width (equal to the 
reduced undulator length), the wave vector fcg and the 
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width of envelope cast into: 

ka - V . 

Au(l + tf 2 )’ 

J_ = a u (i + a^ 2 ) 

CTfc 47<7 e 


where <r e is the dispersion of electron’s energy, deter¬ 
mined in Sect. El 

The reduction of the wave vectors and the envelope 
width, and the displacement r read 


k m = h 0 - mfcgAc/7 ! 


2 2 / 

°m = O’l/W = 


An 


1 5 


to 28fcg yyAc 


, s 27TI/1 , . 

T n -m = (n- to)—— = (n- m)Ti 
ko 


Ultimate coherency - maximum on-axis density of ra¬ 
diation - is attainable by means of keeping the carrier 
frequency constant. This can be done by proper tapering 
of the undulators, see [Iq!, IH, UK- Further we will refer 
to the cases of constant carrier frequency as ‘tapered’. 

Examples of the interference terms are plotted in 
Fig. [6] As it can be seen from the plots, with increase 
in the pulse number to interference with the neighboring 
pulses decreased faster. For a general (untapered) case, 
the interference terms undergone oscillations. 

Energy density T> m as a function of number of emitted 
photons, beginning from the front end to = 1 is plotted 
in Fig. [Tl the coherence factor in Fig. [HI 

As it can be seen from Fig. [7) radiation from the front 
end of the undulator displays coherency. With increasing 
in the number of radiated out pulses radiation becomes 
less coherent. Deviation from the coherent limit (the 
blue curve in Fig. 0 is more severe for the non-tapered 
undulator. 




FIG. 6. Interference terms magnitude / y m ,n for different cen¬ 
tral pulse numbers m = 1,10, 20 ,..., 100, top - ‘tapered’ un¬ 
dulator, bottom - untapered one. 


of the coherent undulator periods that turn out to be 
independent of the undulator parameter K : 


tu-M ^ /~5~ / q A u \ 
coh ~ 47 tV 63 \ jXcJ 


2 x 10 8 


A u [m] 
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(25) 


D. Coherency criterion 


A criterion that indicates coherency length extends of 
the undulator downstream, counted from the frontend 
(the most coherent radiation), may be established from 
a mere contraction that the exponent power in (1221) is 
equal to minus unity: 


71,6. ~ exp(-l) 


cr i cr |» (k\ - %J 2 +• 

2 R+4) 


= 1- (23) 


For the most coherent tapered undulator, the number 
of coherent pulses is 


At) _ t f^( hM 

^ V 63 V?AcJ 


K 2 . 


(24) 


Multiplying £* by number of the undulator periods for 
emitting a single photon iq, Eq. (usd, we get the number 


with A u [m] being scaled in meters. 

The case of untapered undulator is more complicated. 
Nevertheless, for a limiting case when the undulator pa¬ 
rameter is not too small, 


K 2 

1 + K 2 


> 1CT 9 



(26) 


the number of coherent pulses is determined only by the 
undulator parameter: 


^ ut) » 3(1 + A' 2 ) 2 . (27) 


Positions of corresponding coherency criteria are indi¬ 
cated with diamonds in Figs. 17181 As it is anticipated, 
the tapering substantially increases coherency of single¬ 
electron radiation. 



















FIG. 7. Density vs. number of pulses; the blue curve rep¬ 
resents full coherence, the green one is for no coherence, the 
red for a ‘tapered’ undulator, the black one is for a regular 
undulator. The diamonds indicate the coherence criteria. 



FIG. 8. Coherence factor vs. number of pulses, the black 
curve represents a regular undulator, the red one a ‘tapered’. 

IV. SUMMARY AND CONCLUSIONS 

A novel method for evaluating the energy spectrum of 
electrons emitting the undulator- and the inverse Comp¬ 
ton radiation was developed. The method based on Pois¬ 
son weighted superposition of electron states is applica¬ 
ble for whole range of the emission intensity per elec¬ 
tron pass through the driving force, from much less than 
unity emitted photons (Compton sources) to many pho¬ 
tons emitted (undulators). The method allows for ac¬ 
count contributions in the energy spread both from the 
Poisson statistics and diffusion due to recoils. As it was 
shown, for the case of small number of photons emitted 
per pass the electron spectrum is mainly determined by 
the spectrum of recoils, while for many photons emitted 
major contribution comes from the Poisson law of recoils 


distribution. The theoretical results were confirmed by 
simulations. The evaluated width of electron spectrum, 
see o, is larger than that derived in m for ID Poisson 
model or the diffusion model considered in [ 12 ] since it 
accounts for the both processes. 

The electron energy spectrum was used for evaluation 
of the on-axis density of photons and their coherency 
making use of the ‘carrier-envelope’ presentation for the 
emitting photons. The evaluated maximum coherency 
degree of single-electron radiation is obtained to be pro¬ 
portional to the undulator spatial period and inversely to 
the energy of electrons, the number of coherently emit¬ 
ted undulator periods almost independent of the undu¬ 
lator deflection parameter (for the ultimate case of the 
tapered undulators). The results of our study are appli¬ 
cable both for the classical limit of low-energy undulator 
radiation and for the quantum limit of Compton gamma- 
ray sources. 

‘Classical limit’ for the coherency criterion corresponds 
to the case when all of photons are emitted at the reso¬ 
nant frequency and the pulse width (envelop) equal to the 
reduced length of an undulator, cr p = iV u A p (l+A' 2 )/( 47 r). 
The number of photons in the classical limit, h —> 0, 
would be taken equal to the number of undulator peri¬ 
ods. Thus we have the classical density V c \ oc IV 2 , see 

e.g. m- 

The Compton sources produce incoherent radiation. 
This statement follows from a simple estimation: Even 
for a long-period laser, Ai as = 10 //m, and a low-energy 
electrons, 7 = 100 , the number of coherently radiated 
periods would not exceed 10 . 

The introduced coherency criterion C3 does not work 
for number of pulses emitted in average by single electron 
less than unity, £ < 1. Rigorously speaking, due to Pois¬ 
son distribution of photons in the case (C 1 coherency 
modulation in angular distribution will be of order £: 
P 2 (£)/Pi(£) ~ £/ 2 , which is least of practical interest. 

Increase in the undulator brightness, which is directly 
connected to the coherency of single-electron on-axis ra¬ 
diation, for higher energy of emitted photons requires 
enlarging the undulator spatial period together with in¬ 
creasing the electrons energy: A u /q = const. In this 
case the energy of on-axis photons will increase in di¬ 
rect proportion to the energy of electrons. Tapering of 
the undulators are crucial for coherency of radiation, un¬ 
tapered undulators provide much less on-axis density of 
radiation. 
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